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As an extension of the theory of Dyson's Brownian motion models for the standard Gaussian random-matrix 
ensembles, we report a systematic study of hermitian matrix-valued processes and their eigenvalue processes 
associated with the chiral and nonstandard random-matrix ensembles. In addition to the noncolliding 
Brownian motions, we introduce a one-parameter family of temporally homogeneous noncolliding systems of 
the Bessel processes and a two-parameter family of temporally inhomogeneous noncolliding systems of Yor's 
generalized meanders and show that all of the ten classes of eigenvalue statistics in the Altland-Zirnbauer 
classification are realized as particle distributions in the special cases of these diffusion particle systems. 
As a corollary of each equivalence in distribution of a temporally inhomogeneous eigenvalue process and 
a noncolliding diffusion process, a stochastic-calculus proof of a version of the Harish- Chandra (Itzykson- 
Zuber) formula of integral over unitary group is established. 



It is interesting to consider today mathematical-physical sequences of the two classic papers 11 and 10 by 
Dyson of random matrix theory, which appeared sequentially in the same volume of the journal in 1962. In 
one of them 11 , following the early work of Wigner, he gave a logical foundation for his classification scheme 
of random-matrix ensembles based on the group representation theory of Weyl and established the standard 
(Wigner-Dyson) random matrix theory for the three ensembles called the Gaussian unitary, orthogonal, and 
symplectic ensembles (GUE, GOE, and GSE). He introduced in the other paper ^Hl the hermitian matrix- 
valued Brownian motions, which are associated with these Gaussian random-matrix ensembles, and studied 
the stochastic processes of eigenvalues of the matrix-valued processes. Combining the standard perturbation 
theory of the quantum mechanics and a simple but essential consideration of the scaling of Brownian motions, 
he generally proved that the obtained eigenvalue processes are identified with the one-dimensional systems 
of Brownian particles with the repulsive two-body forces proportional to the inverse of distances between 
particles. These processes are now called Dyson's Brownian motion models Y(t) = {Yi{t), Y2{t), • • • , Yjv(O) 
described by the stochastic differential equations 



with (3 = 1,2,4 for GOE, GUE and GSE, respectively, where Bi{t),l < i < N are independent one- 
dimensional standard Brownian motions. Dyson's classification scheme has been extended. In addition 
to the standard three random-matrix ensembles, their chiral versions (chGUE, chGOE, and chGSE) were 
studied in the particle physics of QCD associated with consideration of the gauge groups and quantum 
numbers called flavors 54, 53, 27] H|. After that extension, Altland and Zirnbauer introduced more four 
ensembles called the classes C, CI, D, and Dili for the solid-state physics of mesoscopic systems considering 
the particle-hole symmetry, which plays an important role in the Bogoliubov-de Gennes framework of the 
BCS mean-field theory of superconductivity P l^- These totally ten Gaussian ensembles are systematically 
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argued by Zirnbauer [HEI based on Cartan's classification scheme of symmetric spaces and Efetov's 
supersymmetry theory ^S] ■ 

One consequence of a combination of the two papers by Dyson may be to give a systematic study of 
matrix- valued diffusion processes (i.e. diffusion processes in groups or algebraic spaces) and perform the 
classification of eigenvalue processes as generalization of Dyson's Brownian motion models. This line has 
been taken by Bru [Hid, Grabiner Konig and O'Connell [HSI and others, and one of the purpose of the 
present paper is to clarify the relationship between statistics of (nonstandard) random matrix theory and 
stochastic processes of interacting diffusion particles in the type of Dyson's Brownian motion models studied 
in the probability theory. We will claim in Sec. II that the matrix-valued processes called the Wishart process 
by Bru 7 and the Laguerre process by Konig and O'Connell 36 are the stochastic versions of chGOE and 
chGUE, respectively, in the sense of Dyson and derive in Sec. Ill the diffusion processes describing the 
eigenvalue statistics of the classes C and D of Altland and Zirnbauer, following Bru's matrix- version of the 
stochastic calculus based on the Ito rule for differentials. 

Due to the strong repulsive forces in the processes of the types of Dyson's Brownian motion models, 
particle collisions are suppressed. Impossibility of collision may be generally proved by the same argument 
as Bru, who showed that the collision time between two eigenvalues of the Wishart process is infinite 
(r = -foo a.s.) Jo^. For the (3^2 (GUE) case of Dyson's Brownian motion model (Q}, if Y(0) G 
then Y{t) G for all t > with probability 1, where denotes the Wcyl chamber of type Aat-i; 
= {x e K^; Xi < X2 < ■ ■ ■ < xn}. Using the Karlin-McGregor formula ,2S« JS, the transition density of 
absorbing Brownian motion in W;^ 
the determinant 



the absorbing Brownian motion in W;^ from the state x at time s to the state y at time t{> s) is given by 



f^(t-s,y\x)= det 

^ ^ l<i,j<N 



G^{t-s,yj\x,) 



x,yeW^, (2) 



where each element is the Gaussian heat-kernel G-^{t,y\x) — e ''^ /^*/ \J1-Kt. Grabiner |21j pointed out 
that the transition probability density of the process with /3 = 2 is given by 

P^(s,x;i,y) = _ 5,y|x)/jA(y), 

where h^{x) — rii<i<j<Ar(2^j ~ ^i)- Since h^{x) is a strictly positive harmonic function in W^, this is 
regarded as the ^.-transform in the sense of Doob [HI , and it implies that the eigenvalue process of GUE is 
realized as the noncolliding Brownian motions {i.e. the /i-transform of an absorbing Brownian motion in the 
Weyl chamber of type Ajv_i). Konig and O'Connell also showed that the eigenvalue process of the Laguerre 
process, which corresponds to chGUE, is realized as the noncolliding system of the squared Bessel processes 
. In Sec. IV, we show that the eigenvalue processes of random matrices in the symmetry classes C and D 
of Altland and Zirnbauer are realized as the noncolliding system of the Brownian motions with an absorbing 
wall at the origin j85| {i.e. the /i-transform of an absorbing Brownian motion in the Weyl chamber of type 
Cjv) and as the noncolliding system of the reflecting Brownian motions {i.e. the /i-transform of an absorbing 
Brownian motion in the Weyl chamber of type Djv), respectively. These three kinds of systems are discussed 
as special cases of a family of noncolliding systems of diffusion particles with one parameter > —1, in which 
each particle is following the d = 2{i> + l)-dimensional Bessel process defined by the transition probability 
density [HlUHj 

G^''Ht,y\x) = y^ie-(-W)/2tj^ for x>0,y>0, 

where F denotes the Gamma function and is the modified Bessel function; I,y{z) = J2'^=oi'^ f'^)'^^'^'^ /{^(''^^ 
l)T{u + n+l)}. 

How can we realize other six eigenvalue processes in Altland-Zirnbauer's ten classes of random-matrix 
ensembles as well by noncolliding systems of diffusion processes ? In our previous papers |^ 02] we consid- 
ered the situation that the noncolliding condition is imposed not forever but for a finite time-interval (0, T] 
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to define the temporally inhomogeneous noncoUiding Brownian motions X(i) — (X\[t)^X2{t\ ■ ■ ■ ■,Xj^[ty). 
Of course, we can see that X(t) Y(t) in distribution as T 00. We observed for the finite time-interval 
t e [0, T] that, if we set X(0) = Y(0) = with = (0, 0, • • • , 0) G M^, then 

where C[A] = (27r)^/2 r(i), C[A'] = 2^/^ J],^! r(i/2), and V^[A] = N[N - l)/4. This is regarded 
as a multivariate version of the Imhof relation in the probability theory ,25 , since it implies the absolute 
continuity in distribution of the temporally homogeneous process Y(i) and the inhomogeneous process X(t) 
in [0, r], but from the viewpoint of random matrix theory the important consequence of this equality is the 
fact that the process X(i) exhibits a transition in distribution from the eigenvalue statistics of GUE to that of 
GOE and thus the GOE distribution is realized at the final time t ^ T. In Sec.V, we develop this argument 
by replacing the Brownian motions Xi(t),l < i < N hy the generalized meanders with two parameters 
{v, k), V > —1, K S [0, 2(v + 1)), introduced as the temporally inhomogeneous diffusions associated with the 
Bessel process by Yor I^H] , whose transition probability density is given by 

4^^'') {s, x; t, y) = ^__gM [t ~ s, y\x)hi,^'''^ {t, y) (5) 

for 0<s<t<T, a;,?/>0 with h!^''"\t,x) = dz G^^^T — t,z\x)z~'^. By choosing the two parameters 
(r^, k) appropriately, this family of noncoUiding systems of generalized meanders provides such diffusion 
processes that exhibit the transitions from chGUE to chGOE and from the class C to the class CI. We will 
also consider the processes, in which the noncoUiding condition collapses at the final time t = T in the ways 
that all particles collide simultaneously or only pairwise collisions occur. In the special cases in the latter 
situation, we have the processes showing the transitions from GUE to GSE, from chGUE to chGSE, and 
from the class D to the class DHL 

The present study of the temporally inhomogeneous noncoUiding diffusion processes gives two kinds of 
byproducts, (i) Topology of path-configurations of our processes on the spatio-temporal plane Ex [0, T] is 
determined by the conditions at t = and t = T. We will be able to discuss the topology of random directed 
polymer networks [S] 1141 using the random matrix theory. Such correspondence between the topology of 
path-configurations and random-matrix ensembles is recently used by Sasamoto and Imamura to analyze one- 
dimensional polynuclear growth models 149' . (ii) A variety of versions of Harish-Chandra (Itzykson-Zuber) 
formulae of integrals over unitary groups 22, 26 are derived as corollaries of the equivalence in distribution 
of the eigenvalue processes of matrix- valued processes and noncoUiding diffusion processes. Other remarks 
are given in Sec. VI. 

II BRU'S THEOREM 

A Hermitian matrix-valued stochastic processes 

We denote the space oi N x N hermitian matrices by 7i(iV), the group N x N unitary matrices by U(iV), 
and the group of TV x iV real orthogonal matrices by 0{N). We also use the notations S{N) and A{N) for the 
spaces of N X N real symmetric and real antisymmetric matrices, respectively. We consider complex- valued 
processes £,ij{t) & C,l < i,j < N,t G [0, 00), with the condition ^ji{t)* = ^ij{t), and define the matrix- valued 
processes by E{t) = {^ij{t))i<ij<N G H{N). We denote by U{t) = {uij{t))i<ij<N the family of unitary 
matrices which diagonalize S(t) so that 

U{t)^E{t)U{t) = A{t) = diag{Ai(t), A2 • • • , XN{t)}, 

where {Xi{t)}fLi are eigenvalues of and we assume their increasing order 

Al(t) < A2(i) < ••• < AAr(t). (6) 
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Definery(t),l <i,j < iV, by r„- = (t/(t)tdS(t){/(t))y (t/(t)tdS(i){/(t))j„ where dS(t) = (c^e^J)l<^J<JV• 
We denote by the indicator function: l(w) = 1 if the condition ui is satisfied, and l{uj) — otherwise. 
The following theorem is proved for the stochastic process of eigenvalues X{t) = {Xi{t), \2{t), ■ ■ ■ , ^N{t))- 

Theorem 1 Assume that ^ij(t),l < i < j < N are continuous semimartingales. The process X{t) = 
{Xi{t), \2{t) , ■ ■ ■ , Xn (t)) satisfies the stochastic differential equations 

dX^{t) = dAU{t) +dJ,{t), l<i< N, 

where Mi{t) is the martingale with quadratic variation {Mi)t = Jq Tii{s)ds and Ji{t) is the process with finite 
variation given by 

^ 1 

^^^W = E x^{t)-x,it) ^^^'^*^ ^ x,it))r,,{t)dt + dT,{t) 

where dTi{t) is the finite-variation part of {U{ty d'E.{t)U{t))ii. 



Since this theorem is obtained by simple generalization of Theorem 1 in Bru |H] , we call it Bru's theorem 
here. A key point to derive the theorem is applying the Ito rule for differentiating the product of matrix- 
valued semimartingales: If X and Y are N x N matrices with semimartingale elements, then 

d(xty) ^ {dx^Y + xt(dr) + (dx)t(dr). 
B Four Basic Examples 

Let N = {0, 1, 2, • • •} and assume v £ N. Let Bij{t), Bij{t), 1 < i < N + v,l < j < N he independent 
one-dimensional standard Brownian motions. For 1 < i,j < N we set 



Sy (t) 



1 
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— Bij{t), ifi<j. 



Bii (t), 
1 



if i = j, 
Bji{t), iii>j, 



and Uij (t) 



1 
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—B,j{t), ifi<j, 



0, 



if i = j, 



Here we show four basic examples of hermitian matrix- valued processes and applications of Theorem ^ 

(i) The first example of hermitian matrix-valued process is defined by 

iCij{t))i<ij<N = {sij{t) + \/^aij{t))i<,j<N , t£ [0,oo). 

By definition d^ij{t)d£_ki{t) = SuSjkdt, 1 < i,j,k,£ < N, and thus Tij(t) = 1. Therefore X(t) solves the 
equations of Dyson's Brownian motion model ^ with (3 = 2. 

(ii) The second example is given by 

S(t) = (s„(t))i<,j<Ar G S{N), t e [0,C30). 

In this case d^y(i)d^fc^(t) = {digdjk+SikSje)dt/2, 1 < i,j,k,£ < N, a.ndthusTij{t)dt = {l+Sij)dt/2, 1 < 
hj < N. Then X{t) solves O with /? 1. 



(iii) We consider an {N + v) x N matrix- valued process by M{t) — {Bij{t) 
and define the N x N hermitian matrix-valued process by 



ij )l<i<N+v.l<j<N 



a.{t) = M{t)Hi{t), <e[0,oo). 



(7) 
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Since the matrix :^(t) is positive definite, the eigenvalues are nonnegative. By definition we see that 
the finite-variation part of d£^ij{t) is 2{N + v)5ijdt and d^ij{t)d^ke{t) = 2{^ie(t)Sjk + ^kj{t)Sie)dt, 1 < 
k,£< N, which imply that dTi(t) = 2{N + v)dt and (t) = 2{\i{t) + \j{t)), I <i,j < N. Since 

{Mi)t — /q 4Ai(s)(is, the stochastic differential equations for X{t) are given by 

dK{t)^2^/x;{rjdB,{t)+p\iN + ,y)+ A S-a'S \ ^<^<N, (8) 

with 13^2. 

(iv) Set B{t) = iBij(t))i<i<N+^^i<j<N and define 

= B{tfB{t) G S{N), t e [0, oo). (9) 

We see that the finite- variation part of d^ij{t) is {N + u)5ijdt and d£_ij{t)d£_ke{t) — {£,ik{t)Sji+£,ii{t)Sjk + 
^jkit)Su + iji{t)5,k)dt, 1 < i,j,k,l < N. Then dT,{t) ^ {N + v)dt and T,j{t) = (A,(t) + \j{t)){l + 
Sij), ^ l^i,] ^ N . The equations for \{t) are given by ^ with (3—1. 

The process (jSJ was called the Wishart process and studied as matrix generalization of squared Bessel 
process by Bru 0. Konig and O'Connell 36^ called the process the Laguerre process and studied its 
eigenvalue process (jHl with (3 = 2. 

C Relation with the standard and chiral random matrix theories 

Here we assume that (0) = Bij (0) = for all 1 < i < + i^, 1 < j < iV, and thus the initial distribution 
of 2(t) is the pointmass on an x iV zero matrix O; /j,(S € •;0) = (5o- In this case the distributions of 
S(t)'s are related with those studies in the standard (Wigner-Dyson) random matrix theory j^OI and the 
chiral random matrix theory |54l 1531 ITZl 15 1| . 

(i) Example (i) and GUE. For GUE with variance cr^ = i of random matrices in the space Ti{N) = M'^I^l 
with d[A] = N"^ , the probability density of eigenvalues A in the condition © is given as 00! 

j--d[A]/2 ( |\|2^ 

where |A|2 = Af . For CJ with (3 = 2, pA(0, 0; t, A) = q^^^{\- 1), t > 0. 

(ii) Example (ii) and GOE. The probability density of eigenvalues A with the condition ^ is given as 001 

for GOE with variance ct^ = t in S{N) = d[A'] = N{N + l)/2. If we denote by p^'(s, A; t, A') the 

transition probability density of the process 1^ with /? = 1 from A at time s to A' at time t{> s), then 
p^'{0,O;t,X)=q'^°^{\;t),t>0. 

(iii) Example (Hi) and chiral GUE. We denote by M{N + TV; C) and M{N + v,N;R) the spaces of 
{N + v) X N complex and real matrices, respectively. We see that M{N + N; C) = M2Af(7V+i/) write 
its volume element as V{dM), M e M.{N 4- N; C). The chiral Gaussian unitary ensemble (chGUE) with 
variance t is the ensemble of matrices M E A4{N + v, N; C) with the probability density 

^chGUE(^,. ^ ^^-^^^ exp |--TrMtM| (10) 
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with respect to V(dAI). It is known that any matrix AI e A4{N + h', N; C) has family of pairs (U, V),U E 
U(iV + iy),V e U(iV), which transform M as M = WKV, where K e M{N + v, N; R) is in the form 



K = with X = diag{Ki, K2, • • • 1 «^Af}i > 0, 1 < i < A^, 

and the v x N zero matrix O. We assume that U and V are chosen so that 

< Kl < K2 < • • • < KJV. (11) 

The matrices (J7, if, T^) can be regarded as "polar coordinates" in the space AA{N + v,N\'C). We have 
M^M = yl'Ay, where A = diag{Ai, A2, • ■ ■ , Aat} with the relations Xi = nf,! < i < N. Then k = 
(ki, K2, • ■ ■ 1 ^n) is a set of nonnegative square roots of the eigenvalues of AI^M. Let dii{U, V) be the Haar 
measure of the space \J{N + 1^) x IJ{N) normalized as /u(w+i>)xu(w) ^) = 1 a-^id dK = Hi^i ^'^j- Then 

we can show that 

V{dAI) = i^^^- /x((2-+i)/2) (^fdi^dfiiU, V), (12) 

where = 2^(^+''-i) nti{r(j)r(« + i^)} and 

TV 

l<i<j<N k=l 

For any pair of unitary matrices U G U(iV + 1^) and G U(iV), the probabiHty /xJ;'^°uE(^^^. (-^^j^ 
invariant under the automorphism M ^ U^MV. By integrating over dfj,{U,V), we obtain the probability 
density of k with the condition (jTTll as EHl HIl Elj 

<zf«™(K;t) = L__exp|-^|/z((2''+i)/2)(/.)^ 

Konig and O'Connell studied the process ((SJ with f3 — 2 as a, multivariate version of squared Bessel 
process. Here we consider the multivariate version of Bessel process by extracting the square roots of 
eigenvalues Xi{t) > of S(t) = M{t)^M{t). Setting K,{t) = yj\{t) > 0, 1 < i < in (jHJ with /3 = 2 and 
applying the Ito rule for differentials, we find that n{t) solves the stochastic differential equations 



dZ,{t) = dB,{t) + ^ 



E 



Z,{t) ^\Z,{t)- Zj{t) Z,{t)+Zj{t) 

3 ■Ot^'^ 



dt, l<i<N, (13) 



with {P, 7) — (2, {2i' + l)/2). If we denote the transition probability density of this process by p^'^\s, ■ ■) 
for < s < t < 00, then 



(iv) Example (iv) and chiral GOE. We can see M[N-\-v, N; R) = The chiral Gaussian orthogonal 

ensemble (chGOE) with variance t is the ensemble of matrices B € Ai{N + N; M) C Ai{N + i/, N; C) with 
the probability density 

= (2.).(TV..)/2 exp|-^Tri^^i?| (15) 
with respect to the volume element V'{dB) of A4{N + v, N; R). We can show that 

V{dB) = h^""^ {K)dKdfi{U, V), (16) 
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where d^(U, V) is the normaUzed Haar measure of the space 0{N + iy) x 0(iV) and we have used the notation 
= 2^(^+2--«-i)/27r-^/2n^i{r(V2)r((z + 2iy + l- k)/2)} and thus = 2^(^+'^-2)/2^-w/2 

n,^i{r(V2)r((i + i')/2)}. The probabihty density of k with ^ is given as [51155117711^1 



,-N{N+iy)/2 ( I 12 ^ 



By setting Ki{t) = y/Xi{t), 1 < i < N in ^ with (3 — l,we can show that — (Ki(t), K2(i), • • • , HN{t)) 
solves with (/3,7) = If we denote the transition probabihty density of this process K{t) by 

p^^^'is, ■ ;t, ■) for < s <t < oo, then p^")' (0, 0; k) = q;;'^c;OE^^. ^ > g. 

Ill HERMITIAN MATRIX- VALUED PROCESSES WITH ADDITIONAL 
SYMMETRIES 

A Subspaces of unitary and hermitian matrices 

The Pauh spin matrices are defined as 

l\ f (\ 



which satisfy the algebra = I2, ~ 1, 2, 3, and a^ap = V— 1 ^mpw^w for 1 < /i 7^ p < 3, where /jv 

denotes the N x N unit matrix and e^^pt^ the totally antisymmetric unit tensor. They give the infinitesimal 
generators {X^} of SU(2) by = For > 2, define the 2N x 2N matrices T,^ = In cr^, fi ^ 

1, 2, 3. The matrices {S^} satisfy the same algebra as {cr^}. We will use ao to represent l2- 
We introduce six spaces of matrices as subspaces of TL{2N), 

nf,± i2N) = {He n{2N) : H^ll^ = ±S^i/}, /i - 1, 2, 3. 

It is easy to see that 7i3+(27V) = S{2N) and n?,-{2N) = '/^A{2N). Since we have aheady studied the 
matrix-valued process in S{N) as the example (ii) in Sec. II. B, we will consider here the five subspaces of 
n{2N); ^/^A{2N) and {Hp^(2iV)} with ^ = 1, 2, ct = ±. We also introduce the three subspaces of U(2iV): 

Uo(2iV) = {J7 e U(2A^) : U'^U = Si}, 

U^(2iV) = {(7 e U(2A^) : t/^SpC/ - S^}, ^=1,2. 

The conditions imply that these subspaces, 'Hf_ia{2N) and Up(2A^), have additional symmetries compared 
to 'H{2N) and U(2A^). Concerning the eigenvalues and eigenvectors of the hermitian matrices, the following 
lemma may be easily proved. 

Lemma 2 Assume that fl denotes a diagonal matrix in the form diagjcJi, 012, • • • ,ujn} with uji < 0J2 < ■ ■ ■ < 

(i) Any H G ^/^A{2N) can be diagonalized by U e Uo(2iV) as W HU = 0-3. 

(ii) For ^ = 1, 2 any H E 7i^+(2iV) can be diagonalized by U E U^(2Af) as U^HU = cto- 
(Hi) For n = 1,2 any H e Tlp-{2N) can be diagonalized by U E U^(27V) as HU — VL^a^. 



Remark 

(a) Observing the pairing of eigenvalues in a way, (cji, —uji), 1 < i < N, for \/ ~1A{2N) stated in Lemma 
12 (i), the Gaussian random-matrix ensemble of antisymmetric hermitian matrices was discussed by 
Mehta in Section 3.4 of gOl. 
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(b) The condition for 'U2{2N) addition to the unitarity is equivalent with J = UJU^, where 

J = In^[_^ 

Then \J2{2N) forms the TV-dimensional symplectic group. That is, U2(27V) = Sp(iV,C) n \J{2N). (It 
is called the unitary-symplectic group USp(2iV) in JHI-) The matrices H S ?i2+{2N) are said to 
be self-dual hermitian matrices in the random matrix theory |4U| . The pairwise degeneracy stated in 
Lemma|5|(ii) for 7i2+(2iV) is known as the Kramers doublet in the quantum mechanics. 

(c) The condition for H2~{2N) addition to hermiticity is rewritten as H^J + JH ~ 0, which means 
that H e Ti2-{2N) satisfies the symplectic Lie algebra (see for example ^Hj), that is, H2-{2N) = 
5p{2N, C) n n{2N). Similarly, we can see ni-{2N) = so(27V, C) n H(2iV), where so(2iV, C) denotes 
the orthogonal Lie algebra. We can also see that Ui(2iV) = SO(27V,C) n U(2iV), where S0(2^,C) 
denotes the orthogonal Lie group. 

(d) We can see that n^-{2N) = 7i^_ (2iV), /i = 1,2, where 

ni-{2N) = {^=(4 J^j^tY- H,eH{N),A2eA{N;C)Y 

■H2-{2N) = {^=(4 ^^T^-- Hien{N),A2&S{N-C) 

where S{N] C) and A{N-, C) denote the spaces of the N x N complex symmetric and complex anti- 
symmetric matrices, respectively. Altland and Zirnbauer studied 'H2-{2N) and 'Hi^{2N) as the sets 
of the Hamiltonians in the Bogoliubov-de Gennes formalism for the BCS mean-field theory of super- 
conductivity, where the pairing of positive and negative eigenvalues (w^,— a;i),l < i < N, stated in 
Lemma El (hi) for fi — 1 and 2 represents the particle-hole symmetry in the Bogoliubov-de Gennes 
theory. They called 7^2- (2iV) and Tii-{2N) the sets of hermitian matrices in the symmetry classes C 
and D [SHI [II E], since sp(27V, C) = Cat and so (27V, C) = Dat in Cartan's notations (see [2S|). 

B Representation using Pauli matrices and application of Bru's theorem 

Let B'^j{t), B^^{t), 0<p<3, l<«<j<-^be independent one-dimensional standard Brownian motions 
starting from the origin. Put 

and <^.(t)= V2 (^7) 

[ BP^{t), ifz = j, [ 0, ifi =j, 

with s1^{t) = s^iit) and af^.(t) = -a^^^it) for i > j and define sP{t) = (sf^. (t))i<,.j<Ar e 5(7V), t € [0,oo) and 
a'it) = (af^.(t))i<,,,<jv G A{N),te [0,oo), for < p < 3. 

We can see that the hermitian matrix-valued process given as the first example (i) in Sec. ILB can be 
represented, if we double the size of matrix to 2iV, as E{t) = J2p=o{i^''(^) ^ '^p) + V—^{0'''(t) ® o-p)}. By 
choosing four terms in the eight terms, we define the following four different types of 27V x 27V hermitian 
matrix-valued processes: 

3 

^^.a{t) = Y,{iP,,{t)®ap)en^,„{2N) for p = l,2, a = ±, 

where 

sP{t) if/i=l, P7^3 or A* = 2, p = 0, 



^aP{t) if/i=l, p = 3 or p = 2, p ^ 0, 
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^/^aPit) or ^ = 2, p = 0, 

sP{t) if = 1, P = 3 or = 2, p 7^ 0. 



We apply Theoremnto the five processes, \/—lA{2N) and {^pcr(i)} with fi = l,2,a — zk. The results 
are hsted below. 



(a) y/—lA{2N): Since T.ij{t) = {1 - ((Si)y)2}/2, 1 < i,j < 2N, the equations of nonnegative eigenvalues 
are 

By changing the time unit as i — > 2t, this equation can be identified with 1)13(1 with (/?, 7) = (2, 0). 



(b) Hi+(2iV): Since Tij{t) = {l + , 1 < ij < 2N, the distinct eigenvalues solve Dyson's 
Brownian motion model ^ with /3 = 4. 

(c) T-Li^{2N): We see Tij{t) = |l — ((Ei)y)^} , 1 < j < 2A^. Then the nonnegative eigenvalues solve the 

equations (O with (/?, 7) = (2,0). 

(d) H2+{2N): Since Tij{t) = |l + [[T,2)ij)'^} , f < «, J < 27V, the distinct eigenvalues solve the equations 

^ with /3 = 4 

(e) n2-(2N): We can see ry (t) — |l — ((I]2)ij)'^} , 1 < j,j < 2A'^. Then the nonnegative eigenvalues solve 

the equation lO with (/3,7) = (2, 1). 

C Relation with standard and nonstandard random matrix theories 

(i) The eigenvalues of any matrix in the space 7^2+ (2iV) = E**!^ 1 with d[A"] = N{2N — 1) are pairwise 
degenerated (the Kramers doublets) as A = {uji,uji,uj2,i^2, • • • , wat, cjat). We assume that the N distinct 
eigenvalues are always arranged in the increasing order uji < UJ2 < ■ ■ ■ < ^at- For GSE with variance t, the 
probability density of the N distinct eigenvalues in this ordering is given by 001 

where C[A"] — {2tt)^/^ YiiLi ^(2*)- ff we denote the transition probability density of the process with 
/3 = 4 by p^"(s, • ;i, • ) for < s < t < 00, then p^" (0, 0; i, w) = q'^^'^{u;;t), t > 0. 

(ii) We can see that H2- (2A^) ^ R'^I'^l andHi_(2A^) ^R'^^^^ with (i[C] A^(27V+1) and d[D] = 7V(27V-1). 
The probability densities of the processes 'E.2-{t) and Si_(i) with respect to the volume elements V{dH) of 
n2-{2N) and V'{dH) oini-{2N) are given by 

^-d[C]/2 r 1 1 ^ i-d[D]/2 f 1 

-^;^e.^[-^TrH^j, ,-(i/;i) = exp |--Tr 

where c[C] = 23^/27r^(2A'+i)/2 ^^^^j ^ 2^/27r^(2^-i)/2, respectively. As stated in Lemma El (iii) , the 
eigenvalues are in the form \{t) = {uji{t),—uJi{t), uj2{t), —uj2{t), • • • , uJN{t),—LJN{t)). We will assume that 

0<LUl<UJ2<--- <LON- (18) 

Then we have the expressions for volume elements 

VidH) = ^^h^'ic.fdujdU, V{dH) = Ah^^^fdojdU\ (19) 



9 



where dU and dU' denote the Haar measures of \J2{2N) and Ui(2A^), respectively, normahzed as 
IvA2N)dU = 1 and J^^^,j,^dU' = 1. Here C[C] = - (vr/2)^/2 r(2z) and C[B] = C_,/, = 

(7r/2)^/2n,=ir(2i - 1), and h'^{uj) = h^^\uj) , (uj) = h^^^u;). At each time t > 0, for any U G U2(2iV), 
the probabihty 1/^ {H;t)V{dH) is invariant under the automorphism H WHU for H e 7i2-(2-/V), and for 
any U' € Ui(27V), ^jP {H;t)V {dH) is invariant under the automorphism H HU' for H € ni^{2N). 

Altland and Zirnbauer named these two Gaussian random-matrix ensembles the classes C and D, respectively 
(see Remark (d) in Sec. III. A) The probability densities of the N nonnegative eigenvalues with the 

condition (|18|l are then obtained as 

j--d[lJ]/2 f 1^12 ^ 

g«('^;t)--^exp|-^j/.»(u;f for ft = C, D. 

If we denote the transition probability densities of the processes (I13II with {f3, 7) = (2, 1) and with (/3, 7) = 
(2, 0) by p'"{s, ■ ;t, • ) and p^{s, ■ -jt, ■) for < s < t < 00, respectively, then 

p^{0,O;t,u:) = q\u:;t), t>0 for ft = C, D. (20) 



D Real symmetric matrix-valued processes 

Here after, we denote the hermitian matrix- valued processes and Si_(t) by S'-^(t) and S^(i), respec- 

tively. They are given by 

S^(t) = V^a"{t) (g) ctq + s\t) g) (Ti + s'^{t) g) 0-2 -I- S^{t) g) 0-3, 

S°(i) = V^a°{t) (g o-o + V^a\t) (g CTi + V^a^it) (g cr2 + s^{t) g) CT3. (21) 
Since Up, p = 0, 1, 3, are real matrices and a2 is a pure imaginary matrix, if we define the processes as 

S^'(t) = si(t) (gcTi -f s3(t) g)cr3, S°'(t) = x/^a2(t) (gCT2 + s^(i) f^cTa, (22) 

then E^'{t) G 52-(2iV) and E^' (t) e Si-{2N), where 52- (27V) = {S e S{2N) : S'^S]2 = ^ M'^I^'l 

and Si-(2N) = {S € Si2N) : S^Y,i = -T,iS} ^ M''P'l with d[C'] = A^(A^ + 1) and d[D'] = The 
probability densities of S*-^ (t) and EP (t) are given by 

i-dlC']/2 r 1 ^ f-d[D']/2 f 1 



(^;0 = ^prexp|--Tr^^| , (5;t) = -^^..^--TrS 

with c[C'] = 2^7r^(^+i)/2 and c[D'] 2^/277^'/^ respectively. Set 02(27V) = 0(2Ar) n Sp(2iV;M) and 
Oi(2A^) — 0{2N) n S0(2A^; K) and denote their normalized Haar measures by dV and dV', respectively The 
eigenvalues are in the form X{t) = {u;i{t),—uji{t),uj2it),—^2(t), ■ ■ ■ ,ujpf{t),—ujN{t)). Under the condition 
H18fl . we have the expressions for volume elements 

V{dS)^^^h^{i^)du:dV, V'{dS)^-^^h''{i^)du:dV', (23) 

where C[C'] = Ci/2,1 = UZi^i^) and C[D'] = C_i/2,o = 2(^-2)/2p(^/2) n^7^ r(i). The probability 
densities of the nonnegative eigenvalues with H18|l are given as 

j--d[tJ']/2 ( 1^12 ^ 

/(a;;t) = -^cxp|-^|/.«(a;) for = C, D. 

It is remarked that the random-matrix ensemble with the distributions {S;t), whose nonnegative eigen- 
value distribution is given by i^it), is the random-matrix ensemble in the symmetry class CI studied by 
Ahland and Zirnbauer [TII^IS^. 

By applying Theorem^ we can show that the nonnegative eigenvalues of {t) solve the equations p3l) 
with (/3, 7) — (1, 1) and those of EP (t) the equations ifT^ with (/3,7) = (1,0). If we denote the transition 
probability densities of these processes by p'^ (s, • ; i, • ) and (s, • ; • ) for < s < i < 00, respectively, 
then p«'(0,O;t,w) = /(a;;t), t > for ft = C, D. 
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IV TEMPORALLY HOMOGENEOUS PROCESSES 



Assume that v > -1, and we consider the process Y^'^^it) = {Yf'(t),Y^'''{t), ■ ■ ■ ,Y^'{t)),t^ [0,oo), which 
solves the stochastic differential equations (fKUl with (/3,7) = (2, {2v + l)/2), that is, 



dYl''\t) = dB,{t) 




Yl^\t) 



Yj-'^t) Yr{t)+Yr{t) 



dt, 



(24) 



1 < i < N. Remark that \i v = 1/2 and —1/2, the equation is reduced to 113|l with (/3,7) = (2,1) and 
(/3,7) = (2,0), respectively. The Kolmogorov backward equation (the Fokker-Planck equation) for H24|) is 



where b(x) = (6i(x), 
confirm the following. 

Lemma 3 Set 



^pt'') (s, x; y) = i A^p('') {s, x; t, y) + b(x) • ^^p^'^^ [s, x; t, y), 

,6Ar(x)) with &,(x) = (a/92:i)ln/i((2i^+i)/2)(x). By simple calculati 



ion, we can 



det 

l<i,j<N 



G'^-\t,y,\x,) 



(25) 



Then the transition probability density p^^^^s^yi^t^y) from the state x at time s to the state y at time t{> s) 
of the process {24[ is given by 

p('^)(s,x;t,y) = s,y|x)/iW(y), x, y G 



/i(0)(x) 

Since /i/2(a;) — (e^ — e^'^)/V27ra, I_ii2{x) — [e^ + e^'^)/V27ra', if we set 



4. 



GC(t,y|x) 



/2-Kt 



G^{t,y\x) 



y/2^ 



and /»(t,y|x) = deti<,;j<w[G«(t,2;j| 



= C, D, then 



^(1/2) 



(s,x;t,y) = 



/G(t-5,y|x)feG(3 
/iC(x) 



p(-i/2)(s,x;t,y) = 



/D(t-.,y|x)/zD(y) 
/iD(x) 



(26) 



(27) 



(28) 



The above implies the following. Let = {x e :0<a:;i<X2<---<XAr} and W]^ = {x e 
|xi I < a;2 < • • • < xtv}. The former is the Weyl chamber of type Cjv and the latter of type D^r 18 . Since 
and vanish at the boundaries of the Weyl chambers and W]^, respectively, (128(1 implies that 
the processes Y'^^/^^(t) and Y'^~^/^)(i) can be regarded as the A^-dimensional absorbing Brownian motions 
in and in W^, respectively That is, if Y(i/2)(o) e and Y(-i/2)(o) e W^, then Y^^^/'^\t) e 
and Y'^~^^'^\t) G W]^ for all i > with probability 1. Moreover, we notice that (|27() are the heat-kernels 
of the one-dimensional Brownian motion with an absorbing wall at the origin, and of the one-dimensional 
reflecting Brownian motion, respectively @H1- Then, we can also interpret the process Y'^^/^^(t) as the A^- 
particle system of Brownian motions conditioned never to collide with each other nor with the wall at the 
origin in one-dimension 35 , and the process Y^~^/'^^ {t) as the A'^-particle system of reflecting Brownian 
motions conditioned never to collide with each other. For tt =C and D, define 



AA«(t,x) - / dy/»(t,y|x), x G wj^. 



(29) 

Af^itjX) is the probability that N Brownian motions starting from x G does not collide with each 
other nor with the wall at the origin up to time t, and Af^{t, x) is equal to the probability that N reflecting 
Brownian motions starting from x G W]^ does not collide with each other up to time t, respectively. We 
will show their long-time asymptotics in the next section. We can prove the following, which arc consistent 
with (HH) and 
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Figure 1: Process X(t),t £ [0,T], with X(0) = showing star topology. 



Lemma 4 For ly > —1 with fixed t e (0, oo), assume y G W^. Then 



hm pM(0,x;i,y) ^ ' cxp ^ U»^-+i)/^) (y)^ 



|xHo 



/n particular, if v G N, 



\ 2t 

Imi pM(0,x;i,y)=g^'^GUE(y.^)^ 



one? 



hm p(i/2)(0,x;t,y) = g^^y.^)^ p^-^/^) (0, x; t, y) = g°(y;0. 

|x|^0 |x|^0 



(30) 

(31) 
(32) 



Proof. By definition ((23 with ®, if > 0, 1 < Vi < TV, /('')(t,y|x) = {l/t^)llk=iiVk'^^ /^k) 
g-(|x| +|y| )/2t deti<i.j<7v[^iy(a;i2/j/0]- We can use (jA.2p in Appendix A by changing the variables Xi — > xf/2t 
and j/j ^ yf/^^ evaluate deti<ij<Ar[/^(xiyj/t)] and obtain the asymptotic form of /^"•'(t, y|x), 



^-N{N+2u+l)/2 



n 



l<i<j<Af 
N 



Vi) \VtJ 



(33) 



X n (yi-?/?on2^r^-p{-^|x(i+^(^ 

l<k<i<N m=l ^ ^ ^"^ 

in |x|/\/i 0. Using this form in (|26|l . the limit (|30|1 is proved. | 

V TEMPORALLY INHOMOGENEOUS PROCESSES 
A Star topology 

Using 101 the probability that the Brownian motion started at x e does not hit the boundary of 

up to time t > is given by A/''^(i,x) = J^a dy f^{t,y\x). In the previous papers [2I1E21) we gave the 

asymptotic form 



/^(t,y|x) 



(34) 
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in |x|/Vi ^ and showed that 7V^(t,x) = (C[A']/C[A])/iA(x/Vi) x (1 + 0{\x\/Vt)) as \x\/Vt 0. This 
estimate gives that for x e the noncoUiding probability decays in the power-law as t ^ oo 1211 r?7] ; 
Af^{t,yi) ^ with the exponent tp[A] = N{N - l)/4. (Note that ^ is derived readily by using l|XT|l 

in Appendix A.) For a given T > 0, we defined 

for < s < i < T, x,y e W^. Using (jSH) we showed that as |x| ^ it converges to g^(0, 0;i,y) = 
^'j^m]t-dlM/2/C[A'])e-^y^''/^*h^{y)M^{T - t,y). This function g^(s,x;t,y) can be regarded as the transi- 
tion probability density from the state x e W;^ at time s to the state y e at time t{> s) conditioned to 
stay inside W;^ up to time T and defines a temporally inhomogeneous diffusion process, which we denoted 
by X(t) = {Xi{t), X2{t), . . . ,XN{t)), t G [0,T] in Sec. I. This represents the A^-particle system of Brownian 
motions conditioned not to collide with each other in a finite time-interval (0, T]. The process X(t), t £ [0, T], 
starting from X(0) = is illustrated by Figure 1, whose spatial-temporal path-configuration is said to be in 
star topology in the theory of directed polymer networks [14) . As mentioned in Sec. I, this process exhibits a 
transition of the eigenvalue statistics from GUE to GOE |31l \',V2\ . 

In the present section, we consider the temporally inhomogeneous diffusion process associated with 
Y^'^\t) studied in the previous section. We consider the iV-particle system of generalized meanders (01 
conditioned that they never collide with each other for a time interval [0,T]. The transition probability 
density is given by 



for < s < t < T,x,y e W^, where /-^"''"^ (s, x; y) ^ det i<i,j<N[G'"q^''"\s,Xi]t,yj)] with ^ and 
M!r^^^\t,^) = V dyf^^^^\t,^-T,y). Note that 4'^''') (s, x; y) = /M (t - s, y|x)4^''') (i, y)/^^''') (s, x). 



where /i^^'"^ (i, x) = Jlti h!j:'''\t,x^). Since limt^o C^^^t, z\w) = 6{z^w)l{z > 0), hp''\T,^) = Iljli ^^J'' 
for X G W^, and then H35|l can be written as 



ff(r-«)(.,x;t,y) = -/('')(t-.,y|x)AA(''-«)(r-t,y) (36) 



with 

N 



AA('^'^)(t,x)= / dy/(^Ht,y|x)TTyr. 



(37) 



Lemma 5 Assume that v > —1 and k E [0, 2(j/ + 1)). Let x, y e W^^. 



(i) ForO<s<t<T, lim 5f '''^(s,x;t,y) =p^''^(s,x;t,y). 



^ ..... O^"''^''''' --^ - r.i-')! 

(it) ForO <t <T, 



g^r''^\0,O;t,y) ^ lim g^f'^\o,^;t,y) 



ji-N{N+2v+1-k)/2 



^^v[-''^]h^'''^^\y)U^'^-'^\T^t,y). 

(38) 



rp-N(N+2v+l-K)/2 ( I |2^ 

(m) ForT>0, lim 5^''''^^ (0, 0; t, y) = — exp j /i(2-+i-'^)(y) 



13 




Proof. Using ^ for ijSZJl, we have the estimate of TV^'"^'') (t, x) in |x|/^/^ as 



^-W(Ar+2y+l)/2 



AA(''-«)(t,x) = J] 

i<i<r 

•^''^N l<fc«'<Af m=l 



2 / X 2 
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where we have used a version of Selberg's integral formula |50[ I38| 

/ du n - "^'1'" n i^^f "-^e-i"iv^ ^ 2"^+.A^(^-i) n + - 1)) 

i<i<i<JV fc=l 1=1 ^ 

by setting a + 1 - k/2 and 7=1/2 (see Equation (17.6.6) in gO]). By ^ and (EHll, (i) and (ii) are 
obtained. Since limt^Q G'^'^^t,y\x) — S{y — x)l{y > 0), we have limt_^o /'^^H^i yI^) = Y\^=i^{yi ~ ^i) foi" 
x,y e W^. Then limt_,oA/'(''''^)(t,x) = H^^i a;r''l(x ^ W^) and (iii) is obtained. | 

Now we define the process X('''«)(t) = (X^'''"^ (t), ^^''^''^(t), • • • , X^^'^^t)), t e [0,T], as the temporally 
inhomogcneous diffusion process, whose transition probability density is given by for < s < i < 
T, X, y e and PH| for < t < T, y e W^. This process solves the stochastic differential equations 



dt, tG[0,Tll<i<N, 



where (i, x) = (5/9x,) InTV^"^'") (i, x), 1 < i < iV. 

Here we consider the special cases (i^, k) — (1/2,1) and (i', k) = (—1/2,0). By the definitions (|29|) and 
(EH), A/'(i/2,i)(i^x) =:A/'C(t,x)/n.=ia;* and7V'(-i/2^°)(t,x) =A/'D(i,a;), and then ^ gives 

5^^/''^'(s,x;i,y) = _^c(r^_ , x) /""(^ ' vl^W - t, y). 
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9t 



(-1/2,0) 



for < s < i < r, X, y e W^. That is, we can interpret the process X.^^^'^'^^ (t) as the iV-particle system 
of Brownian motions conditioned never to coUide with each other nor with the wall at the origin in one- 
dimension during the time-interval [0,T], and the process X(~^/^'°)(t) as the TV-particle system of reflecting 
Brownian motions conditioned never to collide with each other during the time-interval [0,r], respectively. 
The asymptotic forms A/''*(t, x) = (C[tJ']/C[tl])/i'*(x/v^) x {1 + 0{\yi\/Vi)) in |x|/Vt-> for (1 =C and D are 
obtained by (|39|l . and thus we can see the power-laws of the noncoUiding probabilities, A/''*(i,x) ^ ^-'/'Wl as 
t^ooiorxe W5v, ft = C and D with the exponents ip[C] = tp[D] = N{N - l)/2. As a corollary of 

Lemma El we have the following. 



Corollary 6 (i) For < t < T , if e 

Ji/2,i)/n n.+ - Z r 



5^^/^'^)(0,0;t,x) = ^^c^exp|-^}/.«(x)AA«(r-t,x), 



("m; ForT> 0, ifxe 



lim g^^/^'^^ (0, 0; x) = (x; T), ^lim g^"^/^'"^ (0, 0; t, x) = (x; T) . 

Figure 2 illustrates the processes X'^^/^'^'(t) and X'^^/^'*'^ (i) both starting from 0. The path-configurations 
are in star topology. In the former any particle can not collide with the wall at the origin, while in the latter 
the leftmost particle is reflected at the wall. Another corollary of Lenima[Slis the following. 

Corollary 7 //i^eN,xeW^, lim 3^"'''+^^ (0, 0; t, x) = q;:'''="="^(x; T). 

t y^T 

The combination of Lemma|Hl(i) with 1)31(1 and 1)32(1 of Lemma^l Corollaries and |3 implies that X(^/^'^)(i), 
X(-i/2,o)(^) and X'^'^-^+'^^t) with i/ e N, aU starting fr om 0, exhibit the transitions from the eigenvalue 
statistics of the class C to the class CI, from the class D to the class associated with studied in Sec. III. D, 
and from chGUE to chGOE, respectively, as time t goes on from to T. (See Theorem |21 below.) 

At the end of this subsection, we discuss the relation between the temporally homogeneous diffusion 
process Y^'^'{t) and the temporally inhomogeneous diffusion process X'^'^'^^t) for t e [0,r]. For a time 
sequence to = < ti < ■ ■ ■ < ti-i < tg = T < 00 with £ G {1, 2, • • •}, we consider the multi-time probabilities 
with the initial state Y('')(0) = X('^''')(0) = x(") 



px(o) |^Y('')(ti) G dx(i), . . • , YM(i,) G dxW j = np('')(t._i,x(-i);i„xW)dxW, 

and 

px(o) I^^M g ^^(1) x(^'-) (t,) e dx(^)^ = f[ 5^^^'') (t.^i, x(^-i) ; U, x('' )dx(*) , 

where we have used the Markov property of the processes. Assume that x'^'^^ = and x^*) e I < i < i- 
We use the formulae (126(1 and 1(36(1 and apply Lemmas 0] and [S] Then we have the equality 

ng(r''=Ht._i,x(-i);t,,xW)=r^(^+«-i)/2-g^npM(t,_i,x(- ^ 



' z— 1 ^ ' 



Since this equality holds for arbitrary time sequence to — < ti < ■ ■ ■ < ti^i < tg = T < 00 with 
£ G {1,2,---}, we can conclude the following. 
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Proposition 8 Assume that i/ > e [0,2(i^ + 1)). // X(''^'')(0) = Y('')(0) = 0, then the distribution of 

the process 'X.^'^''^\t) is absolutely continuous with that of the process Y^'^\t) for t e [0,T] and the Radon- 
Nikodym density is given by 

P{Yi-){-) e dw) " a.«/i("Hw(T)) ■ 

When = 1 and {v,k) — (1/2,1), this proposition gives the Imhof relation between the Brownian 
meander and the three-dimensional Bessel process |25j . The relation stated by I0J [21 and the above 
proposition are regarded as the multivariate generalizations of the Imhof relation. 



B Brownian bridges and temporally inhomogeneous matrix-valued processes 



Assume that e N,0 < T < oo. Let B^^{t), B^^{t), I < i < N + < j < N,0 < p < 3 he independent 
one-dimensional standard Brownian motions. For a given matrix m = {niij + V^l"T-y )i<i<w+i/,i<j<7v with 
mij,fhij g M, let {P^)ij{t : niij), {f3l^)ij{t : fhij),l < i < N + i',! < j < N,0 < p < 3 he the diffusion 
processes, which are the solutions of the following stochastic differential equations: 



{P^),,it:m,,)^B^^{t)- 



T-s 
T-s 



-^ds, 

^-^ds, te[o,r]. 



(40) 



The processes (/3^)y(i : Wy) and {(ix)ij{t '■ rriij) are one-dimensional Brownian bridges of duration T both 
starting from and ending at and m^j, respectively |55| . Next for z'' — {z1,j)i<ij<N G S{N) and 
= {z1^)l<^.3<N G A{N), < p < 3, we set 



(4)..(^:4) 



-^(/3^).,(t: V2, 



(/3; 



and 



-i=(/3^),,(t:V2z|;.), 



0, 



if i: 



if i < j, 



if « = j, 



(41) 



with {s^)ij{t : zf^) = {s^)jt{t : z^J and {a^)ij{t : zf^) = -{a^j,)j^(t : zj'J for i > j, where 1 < ij < N,0 < 
p<3 and t £ [0,T]. We define the matrix-valued processes s^{t : z^) = ((s^).y(i : z|^))i<i.j<Ar G 5(A^) and 
a^(t : P) = ((a^)y(i : 2'f,))i<^.j<jv € ^(^). 

In an earlier paper \33\ . we considered the N x N hermitian matrix- valued process "Bxit) ~ s'^{t) + 
y/^a'i^{t : 0),t e [0,T], where O denotes the N x N zero matrix and s^{t) was defined below lfT7|l . This 
process is the temporally inhomogeneous matrix-valued process realized as an interpolation in duration T of 
the first and second processes given in Sec. II. B. Using the invariance in distribution of the process Sr(t) under 
unitary transformations and our generalized version of the Imhof relation we proved the equivalence in 
distribution of its eigenvalue process and X(t) with X(0) =0. As a corollary of this equivalence, we derived 
the formula for any cr G M, 



dU exp 



U(A') 



2a 



C[A]a' 



d[A] 



det 

h^{x)h^{y) i<i,3<N 



G^it,y,\x,) 



(42) 



where dU denotes the Haar measure of U(A^) normalized as 



U(N) 



dU 



1, Ax = diag{a;i, • • • , xjv} and 



Ay = diagjyi, • • • , ?/7v} with x, y S W;^. This is a stochastic-calculus derivation of the Harish-Chandra 
(Itzykson-Zuber) integral formula |22l l2Hj . In this subsection, we give extensions of this argument. 
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As an interpolation of the Laguerre process and the Wishart process (jSJ , we define the matrix- valued 
process 

SLW(t) ^ AfT(t)^MT(t), <e[0,T], 

where Mrit) = {B°j{t) + y/^{(3°)ij{t : 0))i<i<N+^,i<j<N e M{N + v,N;C),te [0,T], where O denotes 
the (A^ + I') X N zero matrix. Similarly, the interpolations between the processes H21(l and H22|l are defined 

by 

S^(t) = \/^a^(t : O) (8) (To + cri + s'^{t : O) (g) CT2 + s^{t) CTs, 

E^{t) ^ y/^QTit : O) (g) (To + y/^a^it : O) » cri + V^a'^{t) (g CT2 + s^(i) (73, 

in which O denotes the iV x iV zero matrix. Let KF^{t) = (Ki^(t), • • • , n^it)), t e [0, T] be the square 
roots of the eigenvalues of S^W(^) ^-^^-^ Q < {t) < ■ ■ ■ < K^{t) and X\t) = (A}(t), A^(0, • ' ' , >^fit)) be 
the nonnegative eigenvalues of Sy(i) with < A}(t) < • • • < A5Y(i) for jl = C and D. We prove the following 
equivalence in distribution among the temporally inhomogeneous diffusion processes. 

Theorem 9 (i) IfiyeNand X^'^-'^+i) (0) = 0, then K^^^it) ^ X^^'^+^l (t), t e [0,T] in distribution. 

(li) //X(i/24)(o) = X(-i/2,o)(o) = 0, then \^{t) = X(i/24)(t) and X^{t) = X^-^/^-°') (t), t e [0,T] in 
distribution. 

Proof. (i) For a given matrix m = (m^ + V~^f>T'ij)i<i<N+ij,i<j<N , mijjfriij e M, we consider A4{N + 
z/, TV; C)-valued process Mrit : m) = ■ + V-T(/3T)y'(^ ■ fnij))i<i<N+iy,i<j<N ,t e [0,r]. From 

the equations H4()|l . we have the equation 

Mrit : m) = M{t) f ' ~ "^ dg, te[0,T], (43) 

where M(t) = {Bfj{t) + \/—lB^j {t))i<i<]\r+u.i<j<N ■ Let TO[/ and mo be random matrices with distribution 
^chGUEj-..^.) j^j^j ^chGOE(-..2.)^ respectively. Since : C) and 0^)-,j{t : C),t e [0,T] are Brownian 

motions when ^ is a Gaussian random variable with variance T independent of Bf^{t) and Bfj{t), if mu and 
mo are independent of M{t),t G [0,T], 

Mrit : mc/) = M{t), Mrit : mo) = MtC*), * e [0, T] (44) 

in distribution. Moreover, since the distribution of the process M{t) is invariant under any transformation 
M{t) WM{t)V, U G U(7V + v),V G U(A^), the following lemma is obtained by the equation (gS)). 

Lemma 10 For any U e V{N + iy),V e \J{N), WMrit : m)V = Mrit : WmV), t e [0,T] in distribution 

By this lemma, if m and m' in A4{N + v,N;<C) have the same radial coordinates, the processes of radial 
coordinates of Mrit : m) and Mrit : m'),t G [0, T], are identical in distribution. Let S^^(t : m) = M^{t : 
m)MT(t : m). Then the above gives the identification in distribution of the processes of square roots of 
eigenvalues of {t : m) and "E}^ [t : m'),t e [0,T]. Now we denote by /j^(-) the probability distribution 
of the process of square roots of eigenvalues of Elf^ {t : m) conditioned that the square roots of eigenvalues 
of m is K = {ki, ■ ■ ■ , Kpf) with the condition We also denote by P(-) and Pt{') the distributions 

of the processes of square roots of eigenvalues of E(t) = M{tyM{t) and ^^(i) = M^{t)AlT{t),t e [0,T], 
respectively. The equalities H44|l give 

P(.)-/ d«P^(-)gf°"^(/^;T), Pt{-)= f d/«P^(-)gf''°^(«;T). 

Then Pt{-) and P(-) satisfy the same relation as the generalized Imhof relation between 'X.^'^''^'^^^ (t) and 
Y^'^'>(t) obtained from Proposition |H1 by setting € N, k = z/ + 1. Since P(-) is equal to the distribution of 
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the temporally homogeneous diffusion process Y'-''\t) (see (jHU of Lemma lU, we can conclude that Pt{-) is 
identical to the distribution of the process X(''''^+^) (t). 

(ii) The second part can be proved by the same argument as the first part. For given y'',zP G S{N), 
yP, zP e A{N), < p < 3, put F = V^lP «> ctq + «) (Ti + «> 0-2 + «> 0-3 e 7^2- (2iV) and Z = 
(8) o-Q + \/^? ® ai + ^/-Iz^ (g) (72 + (g) 0-3 e Tii-(2N). For these Y and Z, we introduce the 
temporally inhomogcneous matrix-valued processes 

S^(t : Y) = \/^a5^(i : y^) (g) ctq + St(^ : 2/^) «) cri + St(^ : y^) <g cr2 + : y^) «> crs, 

: Z) = ^^a^{t : i°) (g (To + \/^a}p{t : z^) g) di + \/^a^(t : z^) (g) a2 + s^{t : z^) g) 0-3. 

The key lemma Uni of the proof is replaced by the following. 

Lemma 11 For any U G U2(2iV),y e Vii2N), WE^{t : Y)U = E^{t : WYU), and V^E°{t : Z)V = 
E^{t : V'^ZV), t e [0,T] in distribution 

For U =C and D we denote by P^'^ {■) the probability distributions of the processes of nonnegative eigenvalues 
of Sy(t : Z) conditioned that the nonnegative eigenvalues of Z is a; = (wi, • • • , cjjv) with (|18|l . We also denote 
by P^{ ) and P^(-) the distributions of the processes of nonnegative eigenvalues of ^^{t) and t G [0, T], 

respectively. Then we have the expressions, 

P8(.)=/ da;P»^^(.)g«(a;;T), P«,(-) = / da; P»^^(.)/ (u;; T) for ft = C, D. 

Comparing them with the {v,k) — (1/2,1) and {v,k) = (—1/2,0) cases of the generalized Imhof relations 
obtained from Proposition |S1 we have the theorem. | 



As a corollary of Theorem |3 the following integral formulae are derived as proved in Appendix B. 
Corollary 12 (i) Assume G N and x,y G W^. For any a G M. 

/ d/i(C/, V) exp |--i^Tr(if, - U^KyV)\K^ - U^KyV) 



/U(Af+i^)xU(Ar) 



det 



h^") {xjh^") {y) i<i,j<N 



where 



O 



^y-\ o 



with = diag{a;i, a;2, • • • ,X]si}, FCy = diag{?/i , j/2 , • • • ,yAf} O'^d v x N zero matrix O. 
(a) Let '^^ C, D. For x, y G W^ , cr G M, 



U(2Af) 



dU exp <^ -— Tr(Ax - ?7UyJ7)^ 



4ct 



^M--'"" det 
/i!'(x)/i*(y) i<i,i<w 



G»((T2,y,>.) 



where Ax = diag{xi,a;2, . . .,xn} ® (73, Ay = diag{?/i,y2, ■ • ■ g" cr3,U(2iV) = U2(2iV) for ^ ^ C 
and U(2iV) = Ui(2iV) for jj = D. 

They are extensions of the Harish-Chandra (Itzykson-Zuber) formula l|42|l. The formula (i) is found in 
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C Watermelon topology 

Consider the A^-particle system of Brownian motions starting from x € at time t = and arriving at 
z e W;^ at time T > 0, which do not collide with each other during the time interval [0,r]. We denote by 
g^''^{0, x; t, y; T, z) the probability density of the state y at time t G [0, T]. It is given by 



g--(O,x;t,y;r,z) = /^^MW0rMbO^ y G , < G [0, T]. (45) 

/A(r,z|x) 



By using (Pll, we can obtain the limit g^-'^fO, 0; i, y; T, 0) lim g^-'^ro, x; y; T, z). Let cttW = 

|xHO,|zHO 

Proposition 13 For yeW^, g^^'"{0,O;t,y;T,O) = q^^''{y;aT{tf),te [0,T]. 

We denote by X'^'^(i), t e [0, T], the temporally inhomogeneous diffusion process, whose probability density 
is given by the above. Its path-configuration on the spatio-temporal plane is illustrated by Figure 3. Such 
a pattern is called watermelon topology in the polymer network theory |14|. 
For 1/ > — 1, similarly to (|45|l we put 

,('')-(0,x;,y;T,z) ^ /^-Ht,y|x)/(-)(T t,z|y) 

y V , , , / /M(r,z|x) 

for X, y, z e W*^, t e [0, T]. By (jSHJ we have the following x ^ hmit. 
Proposition 14 For > -1, x e W^, t G [0, T], 



2N(N+iy) 



5('')-*(0,0;i,x;T,0) = /i^^^+D/^) (x)^ 



exp 



In particular, if v € N, g('')'^(0, 0; x; T, 0) = qJ;'>C!UE(^^ ^^(^)2)^ 5(1/2).* (0, 0; x; T, 0) ^ qC(x, ^j,(t)2)^ 
anrf 5(-i/2).*(0,0;t,x;T, 0) = qD(x, (JtW^)- 

We note that this expression may be formally obtained by taking k 2{v + 1) limit of (jSHJ- 
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D Banana topology 

For e > 0, we consider a subspace of W^^, B^jY(e) = {x = (xi,X2, • • ■ tX2n) G ■ = X2i-i + e,l < 

i < N}. For X e W^^, we will use the notation x°'^'^ = {xi, X3, ■ ■ ■ , X2N-1) and define J^^'^{t,x; e) = 
/BA„(e) ^y°'^'^ /^(i,y|x). We consider the process, whose transition probability density is given by 

A,b. , ^ /^(t-g,y|x)AAA-^(r-t,y;£) . 
5t (s:x;i,y;e) = J\f^.h(T ~ s x-g) ' ^' ^ ^ ^aW' ^ ^ < ^ ^ ^• 

This is the 2A^-particle system of noncoUiding Brownian motions in [0, T] conditioned that the final state at 
time t = T is in B^jY(e). Using ^^,we have 

9t iO,0;t,y;e) = lim g^' (0,x;t,y;£) = - ' 



(e) 

for y e WAy,i e (0,r]. Since linit^o /^(O, y|x) = Vl=iK^^ " 2/^). hmt_o AAA.b(t, x; e) = l(x e BA^(£)), 
and then for y e 

lim,^'^0,O;t,y;s) = Jai ,A(,),-,IV2T ^ »^^(^))- 

As implied in |41| we can take the limit, (7^''^(s, x; i, y) = lim^^o 5'r''^('S, x; y; e), in the above formulae to 
have 

9t (s, x; t, y) = AAA.b(r - ,s, x) ' ^ ^ 

... N{2N+l)/2 -2A^' /,A 



ff^'°(0,O;i,y) = 



^^■''(O, 0; T, y) = 9°^^= (^y^d'i; |) l(y e B^,), (48) 
for X, y e W^jv. < s < t < T, where Af^'^\t, x) = /„a dy f^'^\t, y|x) with 

/^■'^(i,y|x) = ^.^ det_^^^ fG^(i,y,>.) ^G^(i,?/,>,) 



l<i<2A'a<i<Af 



t 



for x e W^jy and y € W^, and B^^ = {x = {xi,X2,- ■ ■ ,X2n) ■ x"^^ e W^,a;2, = X2^-i,l < « < A^}. 
We define the temporally inhomogeneous process 'K.^'^{t)^t e [0,T] starting from or the state in 
and ending at the state in B2Ar as the diffusion process, whose transition probability density is given by 
(|4t)|) - (l48|) . The path-configuration of N particles in this version of noncoUiding Brownian motions on the 
spatio-temporal plane is illustrated by Figure 4, which we would like to call "banana topology" . Important 
point is that at the final time t = T the particle positions are pairwise degenerated and distinct positions 
are identical in distribution with the Kramers doublets of eigenvalues of random matrices in GSE as claimed 
by (EHl). 

Now we consider a 2A^ x 2N hermitian matrix-valued temporally inhomogeneous process defined by 

■~\,{t) = {s°(t) + V^a°,(<:0)}(8)CTo + {4(^:0) + V^ai(t)}(g)cri 

+ {4(^ : O) + V^a^(i)} ®cr2 + {4(^ : O) + \/^a^(i)} ® 0-3, (49) 

where the elements of the N x N matrices {sf^{t : z^), (a^ : z^)}^^-^ are given by grj. By definition, E^{T) 
distributes with the probability density of GSE. Then the same argument as Theorem may prove the 
following. 
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Theorem 15 Let \{t) = (Ai(t), A2(i), • • • , A2Ar(t)) be the eigenvalues of the process \49j) with Xi{t) < X2it) < 
• • • < X2N{t). //X^^'^(O) = 0, then \{t) = X^^'^(i), t e [0,T] in distribution 

As a corollary of this theorem, we will have the following version of Harish-Chandra formula, which is 
found as Equation (3.46) in |41) . 

Corollary 16 Let x = {xi,X2, ■ ■ ■ ,X2n) e W^^, y = (yi,y2, ■ ■ • ,2/Jv) G W^. For any a eR 

wftere Ax = diag{a;i,X2, • • ■ ,X2n}, Ay = diag{2/i,y2, ■ • ■ ,yN} ® ctq, anrf C2Ar[A] = (27r)^n?5i r(«)- 

It is easy to see by the same argument that the transition probability density given below defines the 
temporally inhomogeneous diffusion process 'X.^'^''^^'^{t), t £ [0,T], i/ > —1, k G [0, 2(i/ + 1)), associated with 
X^'^^''), which shows the banana topology: Let 



/('^)^'^(t,y|x) = det 

l<i<2N,l<j<N 



G^''\t,y,\x,) Gi''\t,y,\x,) 



for X e WC^,y e W%, where ci^'^t^ylx) = id/dy)G<^''\t,y\x), and let AA('''«)'b(^^ ^) ^ j^^ jM.b(^^ yj^^ 
nti for ^eWC^. Then 

J^,'.),b. ^. , ^ /(-)ft-s,y|x).A-('--'-).b(r-t,y) 

')N(4N+4iy-l)rTi2N^j--2N(2N+v) 

g(r'-)'^0,O;i,y) = ^ ;i(2.+i)(y)^-|y|V2*_^(.,.),b(y _ y)^ 

5r^)'^0,0;T,y) = ^ ;,((4--2.+3)/4) (yOdd)4^-|y-T/Ti(y ^ bC^)^ (50) 

for x,y e Wgv,0 < s < i < T, where = 2^(2JV+2,.-i) j-j^^ jp(2i)r(2(i + z/))} and B^^ = {x = 
{xi,X2, ■ ■ ■ ,X2n) ■ x°'^'* e W^,a;2i = X2i-i, 1 < i < iV}- We should notice that H5Q() includes the following 
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special cases. 



Here 



j--2N(N+iy) ( \^\2^ ^ 



l<i<i<Af fc=l 

is the probability density of the N distinct square roots n = K2, • • • , i^n) with Ijlll) of the eigenvalues 
of M'^M conditioned that M is a 2N x 2N random matrices in the chiral Gaussian symplectic ensemble 
(chGSE) with variance t jSlESlEaiSI], and 

j--(i[D"]/2 r l/,|2 1 ^ 

.-v;t)=^^cxp|-^l n ^-^fW 



l<i<j<N fc=l 



with d[D"] = 2iV(2iV - 1), C[D"] = = 22^(^~i) Hj^i T{2i)T{2i - 1) is the probability density of the 

nonnegative and distinct eigenvalues lj — {uji,uj2, ■ ■ ■ , i^n) with (|18|l of AN x 4N matrices in the ensemble in 
the symmetry class Dili studied by Altland and Zirnbauer |56[ III 1^ . (Strictly speaking, it is the Dlll-even 
case. The Dlll-odd case is obtained by setting i/ = l/2, K = Oin (|5()|l .'l The above implies that 'X.^'^'^^'^{t) 
with v & N and X.'^~^^^'^^'^{t), both starting from 0, exhibit the transitions from the eigenvalue statistics of 
chGUE to chGSE and from the class D to the class Dili, respectively, as time t goes on from to T. 

A lengthy but explicit expression for the AN x AN hermitian matrix-valued process corresponding to 
X(-i/2,o),b(^) is given as 
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+s'^^{t : O) (g) (0-2 «> <Jp) + V^a^''{t) (g) (0-3 (g) Up) \ 



+ |s°3(0 (g (cto <S> cts) + s"(i) «) (cti g> (J3) 

+^/^af{t : O) (g {a2 <g 0-3) + sf{t : O) g) (0-3 g) cr3)|, 

t e [0,r], where s^'P{t),s'^''{t : O) G S{N) and a''P{t),af^''it : O) G A{N), t G [0,r], are defined similarly 
to H17(l and H41|l . Identification of its eigenvalue process with X(~^/^'''^'^(t) gives the following version of 
Harish-Chandra integral, 



/ dC/exp(--i^Tr(Ax-C/UyC/)^ 

Ju^(AN) I 40-2 



g2Ar[D]a2^(4Ar+l) 
/lD(x)/i(l/4)(y)4 i<,<2Wa<J<A' 



det 

\N,l< 

for any CT e M, X = {xi,X2, ■ ■ ■ ,X2n) e W^^,y = (2/1,2/2, • ■ ■ ,2/w) e W^, where Ax = diag{xi,a;2, • • ■,X2n}'. 
as, Ay = diag{2/i,2/2, ■ ■ ■ ,yN} g> {13 g> do), and C2n[D] = (^/2)^ H^fi r(2z - 1). 
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VI CONCLUDING REMARKS 



In the present paper we showed that the eigenvalue processes of GUE, chGUE, the class C, and the class 
D are realized by the temporally homogeneous noncoUiding diffusion processes and then the temporally 
inhomogeneous noncoUiding diffusion processes were introduced, which exhibit the transitions in distribution 
from the eigenvalue statistics of GUE to GOE, GUE to GSE, chGUE to chGOE, chGUE to chGSE, the class 
C to the class CI, and the class D to the class DHL They are obtained as the special cases of the noncoUiding 
systems of the Brownian motions and those of Yor's generalized meanders. These inhomogeneous processes 
are identified with the eigenvalue processes of the inhomogeneous matrix-valued processes, some of which 
are regarded as the stochastic versions of two-matrix models studied by Pandey and Mehta ^| ^] as 
demonstrated in |^ . We would like to put emphasis on the fact that in order to prove the identification 
we have not used any results by Pandey and Mehta, but used the generalized versions of Imhof relations 
(Q) and Proposition IHJ. Therefore we can give the proof for the Harish-Chandra (Itzykson-Zuber)-type 
integration formulae as corollaries. The present study suggests several open problems. Here we list up some 
of them. 

(i) It does not seem to be possible to realize the eigenvalue processes of the random matrix ensembles 
different from GUE, chGUE, the class C and the class D by any temporally homogeneous noncoUiding 
systems of diffusion particles. Is it possible to realize them as the temporally homogeneous diffusion 
processes with some conditions additional to the simple noncoUiding condition ? 

(ii) Norris, Rogers and Williams 021 studied other matrix-valued process called Dynkin's Brownian 
motion S{t) = G{t)^G{t) with dG{t) = {dB{t))G{t), where d denotes the Stratonovich differential; 
xdy = xdy + dxdy/2 for continuous semimartingales x,y. They showed that the eigenvalues of 

are also noncoUiding systems and derived the stochastic differential equations similar to Q for the 
logarithms of the eigenvalues. As mentioned by Bru (see Remark 2 in 0), G(t) is a matrix-version 
of multiplicative Brownian motion in a sense, while B{t) is the ordinary additional Brownian motion. 
Can we discuss (the logarithms of ) the eigenvalue processes using the random matrix theory and 
noncoUiding diffusion processes as well ? 

(iii) In the non-hermitian random matrix ensembles, eigenvalues are distributed on the complex plane 
|2()l I12j . Is it meaningful to consider the stochastic version of non-hermitian random matrix theory in 
the sense of Dyson 10 ? 

For the temporally inhomogeneous noncoUiding Brownian motions X(t),t £ [0,T] with X(0) = 0, the 
determinantal expressions for the multi-time correlation functions were determined by Nagao and the present 
authors using the self-dual quaternion matrices |44[ I16L 021 and the scaling limits of the infinite particles 
N ^ oo and the infinite time- interval T ^ oo were investigated 03 [SU]. Recently Nagao reported the similar 
calculation on the process, which corresponds to the process X^^^/^'^^ in the present paper 0S|- Calculation 
of the multi-time correlations for the general process X^"^'**^ (t) is now in progress and the study of the infinite 
particle systems will be reported elsewhere |H4| . 
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APPENDICES 



A SCHUR FUNCTION EXPANSIONS OF DETERMINANTS 



Any sequence /i = (/ii,/i2,- 



, ^J■N, 



•) of nonnegative integers in decreasing order /ii > /i2 > • • • > I^-n > 



is called a partition. The non-zero /i^ in fi are called the parts of /Lt and the number of parts is the 
length of ^ denoted by ^(/i). For each partition n with £{ij,) < N, the Schur function defined by Sp(x) = 



deti<i j<Ar(x^^^^ "')/ deti<i j<Ar(x™ gives a symmetric polynomial of order — Mi ^ variables 

xi,X2, ■ ■ ■ , xn G C Note that the denominator is the Vandcrmonde determinant and deti<ij<N{x^~'') — 
(— l)^'^^~^)/^/i^(x) inniEHHl- We can prove the following expansion formulae of the determinants with 
the bases of the Schur functions EI EE] • 



det 

l<i,j<N 



det 

l<i,j<N 



/lA(x)/iA(y) 



^l■.e{^,)<N 



1=1 



where = l/UZi + iV - i + 1) and b^^'^ = l/{n^i r(/i* + N -i + l)r(iy + fM, + N - i + 1)}. Since 
s^(0) = = 0) with = (0, 0, • • • , 0) G , from the above formulae, we have the following asymptotics 



of the determinants. As Ixl 0, 



det 

l<i,j<N 



det 

l<i,j<N 



iVj 



X (l + 0(|x|)), 



N 



,4=1 



nf=ir(j)r(^ + j) 



x(l + 0(|x|)). 



(A.l) 
(A.2) 



B PROOF OF COROLLARY Ha 

By (jSHll of Lemma El(ii), 

gy"^' (0, 0; y) = cxp 



2t 



/i(2^+i)(y) 



dz det 

l<i,i<JV 



r-t'^^'^l 2{T -t) \T -t 



y? + A 



N 

n 

k=l 



rpN(N + K.-l)/2-^-N(N + v) /rp\N(v + l-K) 



X det 

l<i,j<N 



exp 



T 



2t{T - t) 



t{T-ty T 



W / ^ \ V + l-K 



Setting (i/T)zj = a^, 1 < i < TV, t{l - t/T) = cr^ and T/t^ = a, we have 



^~2N ^N(N+2u-K+l)/2 



dae-"l''l /2 det 



N 
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Proof of (i). We write the transition probability density of the process Mxit) by QT{s,mi;t,m2), < s < 
t < T, for mi,m2& M{N + v, N; C). Then by Theorem El (i) and the fact 



(e) \N(N + v) c 

(r'''+i)(0,O;t,y) = l^%^ /,((2.+i)/2)(y)2 / d^,{U,V)QT{Q,0■,t,U^K^V). (B.2) 



9t 



V(N+v)xV(N) 



Weintroducethe7W(7V+i/, AT; C)-valued process M(i)(t) = {rrilht))i<i<N+u,i<j<N and the 7W(iV+j/, iV 



valued process M^'^\t) = {m[f{t))i 



<i<N+i',l<j<N , 



mg)(t) = i?°(t)-|s°(T) + x/^(/30),,(i) and m|f (t) = (T). 

Then Mrit) = M^^^t) + M'^'^\t). Note that {B.° (t) - (T)} arc Brownian bridges of duration T 

starting at and ending at 0, which are independent of {t/T)B°^{T). Hence M'-^^t) is in the chiral GUE 
distribution and M(^^(t) in the chiral GOE distribution, where M^^^t) and M^^''{t) are independent from 
each other. Since E[m'i\t)^] = and E[mlf{t)^] = Qt(0, O; t, M) for M e M{N + v, N; C) can be 
written as 



whose elements are defined by 



I M{N+iy,N;R) 
^N{N+u)/2^-N{N+u) 



„,iV(iv+i/)/2^-iV(iv+iy) r 

, , / rfa;jM(a)e-"l«l'/2-Tr(Af-/f.)t(Af-/^.)/2.^ 



(B.3) 



where we have used the fact H16|l and the formulae H1Q|I . H15|l . Combining ljB.l|l with k = u + 1, l)B.2|l and 
(|B.3|I . we have 



exp ■ 



yf + a] 



U(A'+i/)xU(A') 



Since, for each cr e R, this equality holds for any a > 0, we have the formula (i). 

Proof of (a). By setting [v^k) = (1/2,1) and (z^, k) — (—1/2,0) in (|B.1(I we have the expressions for 



g^^/^-^)(0,O;t,x)= fc^(x) / dag-^l-l'/^ det 

/2 /• 

5r^/^'°'(0,O;t,y) = ^/.-(y) / dae^l^lV^ det 



Following the same argument with the proof of (i) and using the equalities H19II and 123|l . the formulae (ii) 
are proved. | 
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